Coiling of Cylindrical Membrane Stacks with Anchored Polymers 
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We study experimentally a coiling instability of cylindrical multilamellar stacks of phospholipid 
membranes, induced by polymers with hydrophobic anchors grafted along their hydrophilic back- 
bone. We interpret our experimental results in terms of a model in which local membrane curvature 
and polymer concentration are coupled. The model predicts the occurrence of maximally tight 
coils above a threshold polymer concentration. Indeed, only maximally tight coils are observed 
experimentally. Our system is unique in that coils form in the absence of twist and adhesion. 



PACS numbers: 87.16.Dg, 68.10.-m 

The coil motif is ubiquitous in a wide range of natural 
contexts. One-dimensional filaments of mutant bacteria 
101 , supercoiled DNA molecules , and tendrils of climb- 
ing plants |3j all exhibit a writhing instability as a result 
of forcing or interaction with an external agent. Such sys- 
tems are dominated by elastic properties and the appear- 
ance of coils is a result of the relief of twist. In this paper 
we show that coiling can also be effected in cylindrical 
multilamellar tubes of phospholipid bilayers, by anchor- 
ing hydrophilic polymers with hydrophobic side groups 
grafted along the backbone. This system is unique in 
that, in contrast with the above examples, fluid mem- 
branes cannot support any twist. Yet coils are formed in 
the system, and are stable for a very long time. 

Our system is representative of a wide class of systems 
of membranes with embedded inclusions, such as biolog- 
ical cells with membrane-associated proteins. Other ex- 
amples include erythrocyte ghosts incorporating amphi- 
pathic drugs Q and liposomes with covalently attached 
polymers used for drug delivery (for additional ex- 
amples see Q). Our experiments give us a unique op- 
portunity to study a wide variety of phenomena induced 
by such inclusions in a relatively simple and controlled 
environment. 

The cylindrical multilamellar tubes, called myelin fig- 
ures, consist of a large number of bilayers reaching al- 
most to the core [^0- Adjacent bilayers are separated 
by thin hydration layers. Coiling of myelin figures has 
already been observed in phospholipid binary mixtures 
PI when the constituents undergo a phase separation 
process triggered by the addition of Ca++ . Coiled myelin 
figures of egg- yolk phosphatidylcholine (egg-pc) have also 
been reported |l^,|l^. In both cases |]llHl3| coiling was 
attributed to surface adhesion. The novelty of our work 
is that coiling occurs in the absence of both adhesion and 
twist. Our experiments clearly show that surface adhe- 



sion is negligible in our system. 




FIG. 1. A single-stranded coil. This structure formed one 
hour after hydration in a solution containing 0.02 mg/ml poly- 
mer, at a rate of 0.15 pitch/min. Phase contrast images a slice 
of the coil, showing bilayers filling up almost the whole tube, 
though a narrow water-filled core can still be observed along 
the tube axis. The bar represents 10 ^m. 

We believe that in our system coiling results from a 
coupling between the polymer concentration and local 
membrane curvature, induced by anchoring of the poly- 
mer in the membrane. We show below that a high enough 
polymer density together with the constraints imposed 
by the cylindrical geometry of the tube destabilizes the 
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straight tube. Coiling as a manifestation of this instabil- 
ity is proposed here for the first time. 

The coupling between membrane curvature and poly- 
mers has been considered both theoretically JT^-llTt and 
experimentally |]7|,p|,p^-pl| . The theoretical studies em- 
phasize the polymer backbone and its effect on the elastic 
properties of the membrane. They do not consider the 
mobility of anchors embedded in the bilayers. We be- 
lieve this effect is important in our system, and attempt 
to capture it in the model presented below. 
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FIG. 2. Formation of a double helix. After bending, the 
tip slides along the tube, coiling at the same time. The four 
frames were taken at times (a) 0, (b) 20, (c) 30 and (d) 40 
minutes after the double helix was first seen. The bar repre- 
sents 20 /im. 



In our experiments, tubular membrane stacks were 
made of stearoyl-oleoyl-phosphatidylcholine (SOPC), 
with Ci8 alkyl chains. The polymer we used is hy- 
drophilic dextran (MW 162,000 g/mol) functionalized 
both with Ci6 alkyl chains and dodccanoic NBD chains 
as fluorescent markers. The hydrophobic anchors, dis- 
tributed statistically along the backbone (about 1 alkyl 
chains for 25 glucose units) are Cig long. On average 
there are 4 persistence lengths between consecutive an- 
chors. Therefore, the extension of each polymer molecule 
on the two-dimensional membrane is much larger than its 
extension into the third dimension. Samples were pre- 
pared by drying a 0.5-1.0/il droplet of SOPC dissolved 
in a 4:1 chloroform-methanol solution (7.35mg/ml) on a 
glass slide. The sample was then closed, and hydration 
was effected by injecting a polymer solution of known 
concentration, Cp, into the cell. The development of 
myelin structures and their coiling were followed using 
phase contrast microscopy and recorded on video. Our 
experiments were conducted at room temperature, well 
above the solid-liquid transition for SOPC. 

For small values of Cp we observe myelin figures, which 
display a clear tendency to straighten over lengths many 
times larger than their diameter. As Cp is increased, 
myelin figures become more floppy and curved. For large 
enough values of Cp, a writhing instability sets in and 



tubes bend, forming irregular structures, single coils (Fig. 
1) and double helices (Fig. 2). We emphasize that all the 
coiled structures we observe are maximally tight as they 
form and do not tighten up gradually (see Fig. 2); no 
loose coils have been found (unhke Sakurai et al. p^ , p^ ). 
In quantitative terms, this means that the curvature of 
the tube central line, C, is C ~ where rg is the radius 
of the tube. Fig. 2 also demonstrates that there is no ad- 
hesion between membrane surfaces in our system; tube 
segments that are in contact in Fig. 2b are separated in 
Fig. 2c. 

We stress that while the polymer concentration in so- 
lution, Cp, is known, we do not control the surface con- 
centration on the bilayers. The slow evolution of some 
of the structures we observe is consistent with a possible 
variation of this concentration over time. As a first step 
towards a theoretical understanding of this system, we 
neglect this slow evolution. 

We now make the following assumptions: (i) anchors 
penetrate the membrane to a depth of about half a bi- 
layer. This is because the anchor length is compara- 
ble to that of a lipid, and the large hydrophilic back- 
bone to which anchors are attached cannot penetrate a 
bilayer. We assume that these anchors and the poly- 
mer backbone induce a local spontaneous curvature, -ffg 
P,p|j2^ . [it) Polymer molecules are present everywhere in 
the system, including inner regions of the myelin figures. 
We have carried out fluorescence microscopy experiments 
and found that the polymer is present between layers of 
the stack in significant concentrations (additional details 
of the experiments will be published elsewhere), (iii) 
Given that bilayers are in a liquid-like state, polymer 
molecules can diffuse on the membrane. Thus, the en- 
ergy of the system is lowered if they migrate to regions 
where the mean surface curvature H is closer to Hq. (iv) 
Bilayers in a myelin figure maintain a constant area and 
are strongly constrained by the geometry of the stack, 
which prevents them from buckling. We therefore regard 
the myelin figures as flexible cylinders having a fixed cir- 
cular cross section everywhere along their axis. These 
assumptions are strongly supported by our experimental 
observations. Note that the central line of the tube can 
bend. However, as a consequence of the above assump- 
tions, its length remains fixed. 

Based on these assumptions, we present a simple model 
which captures the essential physical features of the ex- 
perimental system, and accounts for many of the exper- 
imental observations. We represent each bilayer as two 
square lattices (in the spirit of lattice-gas models), cor- 
responding to the outer and inner monolayers. We asso- 
ciate two degrees of freedom with each lattice site, which 
corresponds to a membrane patch of area c? . The first 
is the local mean surface curvature, H . The second is a 
binary occupation variable, which takes the values one 
or zero when the site is occupied or unoccupied respec- 
tively by an anchored molecule. This molecule induces 
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a local spontaneous curvature, Hq- The area of a site, 
a^, and the value of Hq depend on the specific mecha- 
nism responsible for the spontaneous curvature. Thus, if 
Hq is induced by individual anchors, the area of a site 
is microscopic (a^ ^ 60A ). If, on the other hand, the 
spontaneous curvature is induced by the polymer back- 
bone, the area of a site is mesoscopic (a^ ^ 6 • 10^ A ). 

By convention, the curvatures of the inner and outer 
layers have opposite signs at the same position. Within 
the model, the energy of the system is a sum of the cur- 
vature energies of the individual area patches: 2kH^ for 
a vacant site, and 2k [H — Hof' for an occupied site. 
K and k' are the local bending rigidities of a single layer 
without and with an anchored molecule, respectively. We 
suppose /v > k; this is consistent with models of com- 
posite membranes |l^Jl^,^ (although the systems these 
models describe arc different from ours). 

In order to find the equilibrium state of a tube, we have 
to calculate its free energy. This free energy depends on 
the curvature of its central line, C, and on p, the average 
of the occupation variable (pa^ is the average density of 
the anchored molecules). If the spontaneous curvature, 
-ffo, is large enough, the free energy of a bent tube is lower 
than that of a straight one. To show this we evaluate 
separately the energy and the entropy of the system. 

Consider one cylindrical bilayer of length / and circu- 
lar cross section of radius r, with the same average occu- 
pancy, /9, on both sides. Let us calculate the energy cost 
of bending the bilayer into a portion of a coil with cen- 
tral line curvature C in two steps. First, the energy of a 
bent cylindrical bilayer with a homogeneous distribution 
of anchored molecules is Ehom — E^om + ^Torm where 
-^hom^ are the energies of the outer and inner monolay- 
ers. According to our model 

Er:n.C)^2pK' j dA[H{C)-Hof 



+ 2(1 -p)k / dA[H{C)] 



(1) 



where H{C) is the local membrane curvature, and is 
known for a cylindrical geometry. The expressions for 
^hom ^hom different since the curvatures of the 
inner and outer monolayers have opposite signs. Note 
that when C 7^ 0, H{C) varies around the bent cylinder. 
For our geometry the total mean curvature J dAH(C) = 
2t:1 for the outer monolayer, while for the inner one 
J dAH{C) = —2111, independent of the central line cur- 
vature C . Thus the cost of bending the cylindrical mem- 
brane, keeping the distribution of anchored molecules ho- 
mogeneous is AEhom — 4[/9k' + (1 - p)k] J dA [iJ(C)]^, 
independent of the value of Hq. Note that the term con- 
taining Hq is independent of C and therefore does not 
contribute to the cost of bending. 

Next, we take into account inhomogeneities in the dis- 
tribution of anchored molecules around the tube. Such 
inhomogeneities reduce the energy if these molecules 



move to regions of membrane curvature closer to Hq in 
both the outer and inner monolayers. The full calculation 
shows that the energy gain, AEinhom{C , Hq) , depends 
linearly on the spontaneous curvature and can become 
arbitrarily large for large values of Hq . 

As for the entropy of the system, we assume that the 
dominant contribution is the entropy of mixing of oc- 
cupied and vacant sites. This entropy is larger when the 
distribution of anchored molecules around the cylindrical 
bilayer is homogeneous, favoring a straight tube. How- 
ever, it does not depend on the spontaneous curvature. 
Therefore, if Hq is large enough, the energy gain due to 
AEinhom{C, Hq) is larger than the free energy cost com- 
ing from AEhom and the entropy of mixing. In this case, 
the tube is bent at equilibrium. It remains to be shown 
that such an equilibrium state can occur for reasonable 
and physical values of the model parameters. 

We have carried out the full calculation of the free 
energy as a function of the central line curvature, the av- 
erage occupancy and the spontaneous curvature (details 
will be published elsewhere) . In the limit of a thick tube 
(a realistic case), the calculation can be done analytically, 
and the free energy of the entire tube (summing over all 
the bilayers) reads: 



F{C,p) 
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where Ktube (p) is the effective bending rigidity of the en- 
tire tube, and I is its length. 

Fig. 3 shows the typical dependence of Ktube on p for 
large enough values of Hq (solid line). When p < p-, 
Ktube is positive and decreases with p. In this regime 
the minimum of the free energy is at C = 0. Therefore, 
the tube is predicted to be straight on the average, but 
with enhanced fluctuations due to the smaller bending 
modulus. Although anchored molecules increase the local 
bending rigidity of the membrane (k' > k) , their mobility 
makes it easier to bend the tube. 

For p_ < p < p+, Ktube is negative, and the free energy 
decreases upon bending with its minimum at the maxi- 
mally possible central line curvature C = 1 /ro (in agree- 
ment with the qualitative argument outlined above). 
Hence, p_ is a threshold occupancy above which straight 
tubes are unstable and form maximally tight coiled struc- 
tures. Above p+ straight tubes become stable again; 
however, this regime is probably unreachable in our ex- 
periments, since too large a polymer concentration de- 
stroys the bilayers. 

A salient feature of our model is that it precludes 
loose coiled structures; i.e. tubes are either uncoiled 
(p < p- and p > /0+) or maximally tight and coiled 
(p_ < p < p+). This is consistent with our experimental 
observations. In addition, according to the model all the 
coils with C = I/tq are equally probable. One should 
therefore expect to see irregular coils without a well de- 
fined chirality as well as regular helices, all of which we 
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indeed observe. 
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FIG. 3. The effective bending modulus of the tube, ntube, 
is parabolic in the average occupancy, p. We have used the 
following values of the parameters: k = lOfcsT, k — 2k 
and ro = 5^m. We find that Ktube depends on a and Hq 
only through the product aHo- The solid curve represents 
Ktube for aHo = 0.3. The physical meaning of this number 
depends on the specific mechanism responsible for the spon- 
taneous curvature. For example, if the spontaneous curvature 
is induced by the polymer backbone, a ^ 8OOA. For the low 
polymer concentration in our experiments (i.e., in the vicinity 
of p = p_), this corresponds to a membrane spontaneous cur- 
vature pHo ~ 0.3p~^ . When aHo is large enough {aHo > 0.19 
for the values of k and k' we have used) Ktube < between 
p_ and p+. For smaller values of aHo, Ktube > for all values 
of p. The dashed curve corresponds to aHo = 0.16. 

A somewhat similar instability of a flat membrane due 
to coupling between membrane shape and local sponta- 
neous curvature has been discussed by Leibler p3| and by 
Safran ||2^. There is, however, a fundamental difference 
between these models and ours. They considered an un- 
constrained flat membrane, whereas our membranes are 
severely constrained by the cylindrical geometry of the 
tube. These geometrical constraints play a crucial role 
in determining the shapes of the observed structures, and 
the final state of the system. 

We believe coiling occurs mainly due to the sponta- 
neous curvature induced by the anchored molecules as 
well as their mobility. Our model emphasizes these as- 
pects, and neglects others such as interactions between 
polymer molecules. We do not expect these effects to 
change the qualitative behavior of the system. So far, we 
have not been able to identify the precise mechanism by 



which the polymer induces spontaneous curvature. We 
intend to study polymer length effects on the coiling phe- 
nomenon in order to partially address this issue. 
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